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Abstract -The Casimir force between graphene sheets is investigated with emphasis on the efTect 
from spatial dispersion using a combination of factors, such as a nonzero chemical potential and an 
induced energy gap. We distinguish between two regimes for the interaction - T — ^ K and T 7^ 
K. It is found that the quantum mechanical interaction (T = K) retains its distance dependence 
regardless of the inclusion of dispersion. The spatial dispersion from the finite temperature Casimir 
force is found to contribute for the most part from n = Matsubara term. These effects become 
important as graphene is tailored to become a poor conductor by inducing a band gap. 
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Introduction . — Long-ranged dispersive interac- 
tions, such as the Casimir 1 force or the van der Waals 
[2] force, originate from electromagnetic field fluctuations. 
They are present between all types of objects at any dis- 
tance, regardless of whether there are permanent electric 
and/or magnetic moments. Despite their universal na- 
ture, the dimensionality, dielectric and magnetic response 
properties, and geometry of the interacting systems can 
influence the dispersive force in profound ways, thus offer- 
ing possibilities to tailor the sign, magnitude, and distance 
dependences of the interaction [3j[4]. 

Single layers of graphite have been recently isolated 
[5j[6], and their dispersive forces present not only chal- 
lenging theoretical problems, but they are also impor- 
tant for graphene based devices [7j , ^ . The electro- 
magnetic fluctuation forces in graphene have been con- 
sidered both at low and high temperature regimes as well 
as with the inclusion of a chemical potential or a band 
gap [£-17 . A distance power law has been obtained 
for the graphene/graphene Casimir force at the quantum 
limit, where the force per unit area(F) goes as the inverse 
distance((i) to the fourth power, F ^ d~^. At the same 
time, due to the graphene gapless electronic structure, the 
fluctuation forces are thermal in nature at room temper- 
ature. This is in stark contrast to the Casimir forces in 
most matter whereupon quantum fluctuation effects are 
dominant even at relatively high temperatures. 



Concerning Casimir interactions in general, the inclu- 
sion of spatial dispersion in the response properties of 
the materials has received significant attention recently, 
since it was shown that the dispersion may play an impor- 
tant role 18j[19j. In this work we investigate to what ex- 
tent spatial dispersion, in combination with modifications 
in the graphene electronic structure through a nonzero 
chemical potential or an induced energy gap, affect the 
graphene-graphene Casimir interaction. In practice, a fi- 
nite chemical potential is induced by doping [2Q{[2T] or by 
the application of external fields 5 . Energy band gaps 
on the other hand can be brought about by the growth of 
graphene on certain substrates [22 , via adsorption 23 or 



hydrogenation 24 . The response properties of graphene 
for the calculations here are taken into account via the 
conductivity obtained from the linear response to the elec- 
tromagnetic field [25^ , from which both the small and large 
spatial dispersion forces are investigated. 

Graphene Conductivity. — An important factor in 
the Casimir force is the dielectric response properties of 
the materials involved. For the graphene/graphene sys- 
tem, these properties are described via the optical con- 
ductivity calculated within the Kubo 26 formalism. Con- 
sider the case when there is a band gap A induced at the 
graphene Fermi level. Near the K-point of the Brillouin 
zone, assuming small wave vector excitations, the Hamil- 
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tonian can be written as 



H 



(1) 



where p is the two-dimensional momentum operator and 
'^0 = 37oa/2 ^ lO^m/s {a = 0.142 nm is the CC distance 
and 7o = 2.4 eV is the nearest neighbor tight binding 
overlap integral [27{[28] ). The Hamiltonian, Eq.Q, corre- 
sponds to a two-band model with eigenvalues and eigen- 
vectors given respectively as: 
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y are unit 



where the wave function is '^(r, s) = |k, 
dimenional wave vector is = -\- kyf {x 
vectors), e*^ = {kx + iky)/k^ Eq = vohk, and s = ±. If 
A = 0, one recovers the gapless Dirac-like model suitable 
for perfect graphene. Using this two-band model and ap- 
plying linear response theory, the dyadic two dimensional 
conductivity is found [25] to be 



see 



eq. Q 

where fo{Es) = l/(e'^^^~^^/^^^ - 1) is the equilibrium 
Fermi distribution function, /c^ is the Boltzmann constant 
and T is the temperature. The velocity matrix elements 
are Vss' = (-5, k|v|k, s'), where v = The first term in 
Eq.Q accounts for intraband transitions, while the second 
one accounts for interband transitions. From hence forth, 
imaginary frequency will be used, which is the relevant 
quantity for calculating Casimir forces. 

In the limit when hvoq 0, the optical conductivity of 
graphene, ao{iuj), is isotropic and it is found as: 
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where /3 = l/ksT. One notes that when fi = and A = 0, 
for ksT <^ hu^ ao{iuj) acquires a universal value ao{iuj) ^ 



(Jo = e /{Ah). In Fig(l(a)), results are shown for ao{i(jj) 



as a function of frequency. The largest differences appear 
at small cj, where the conductivity is most significantly 
reduced when A 7^ (/i = 0) and augmented when /i ^ 
(A = 0). For larger cj, the conductivity approaches ctq. 



A=0.05 eV,n=0 eV 
A=0 eV,n=0.05 eV 
A=0.05 eV,n=0.05 eV 
A=0 eV,n=0 eV 





Fig. 1: (a) The isotropic conductivity ao{iuj) as a function of 
frequency (in eV) for different values of A and /x at T = 300 
K. (b) cTcccc (icc;) , as a function of q for finite A and /x with 
huu = 0.162 eV corresponding to T = 300 K. 



If spatial dispersion is taken into account, ^ is not 
isotropic. For graphene, however, we find that axx and 
(jyy have very similar magnitude and behavior over large 
wave vector and frequency regions, therefore only cTxx is 
shown in Fig( |l(b)| ). The conductivity, axx^ starts changing 
in a more pronounced way for q > 10^ cm~^ . It is noted 
that axx is the conductivity in the direction of the electric 
field parallel to the plane composing graphene sheet while 
(Jyy is the direction perpendicular to the electric field. 

Graphene Casimir Force and Dispersion. — The 

Casimir force per unit area |3 between two graphene 
sheets separated by a distance d is 

see eq. (|8| 

where h{i\u;n\) = iy^e{i\uJn\)lJ^{i\uJn\){^n/c)'^ + q'^ with 
Un = 27rnkBT/h and q is the two-dimensional wave vec- 
tor in the x?/-plane. pe,b are the generalized reflection 
coefficients due to the transverse electric (E) and mag- 
netic (B) field modes, which are given in terms of the two 
dimensional conductivity by 



27T(jyy{iuj, q)uj/[h{iuj)c^] 
+ 27rayy{iuj^ q)u/[h{iuj)c^ 
27r<Jajaj(icj, q)h{iuj)/uj 
. + 27raxx{^^^ q)h{i(jo)/(jO 



(10) 



Thermal Limit. The thermal fluctuation forces are 
first considered, which correspond to the zero Matsub- 
ara term (n = 0) of Eq.([8| in the strong dispersion limit. 
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V(cj,q) = X 



v,,v,, (/o[^.(k - q/2)] - /o[.E,(k + q/2)]) 
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where ficj <C /c^T, fi'L'o<?- In most materials purely thermal 
fluctuation effects are important at distances larger than 
their characteristic thermal wavelength, At = hc/{kBT) 
7 /im. Instead, due to the Dirac-like Hamiltonian of 
graphene and its two dimensional nature, the characteris- 
tic thermal wavelength is significantly reduced - ~ At/200. 
As a result, thermal fluctuation forces become important 
at relatively small scales 16 , d ^ 30 nm. Therefore the 
quasi-static response is a reasonable approximation for a 
large distance range. 

In the quasi-static response limit of hu <^ hvoq <C ksT^ 
the a^x is 

cr^^(za;,q) = (11) 

This corresponds to the longitudinal plasmon excitations 
in the graphene static dielectric function e{q) = 1 -\- Qs/q 
obtained via the Random Phase Approximation(RPA) 
method |29 . For gapped graphene 12 , 
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where a ~ 1/137 is the fine structure constant. Using 
Eq.(|8| and Eq.(ll) the thermal fluctuation stress between 
two graphene sheets is obtained to be 
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One finds that as qsd ^ 1, the stress approaches the 
thermal stress between two perfect metals, i.e.. Ft = 
^BTC(3)/(8^d3). 

A similar calculation can be executed for the problem 
of finding the thermal fluctuation force f{d) between a 
graphene sheet and a polarizable atom 30 with polariz- 
ability x(^). 



see eq. (14) 



Using Eq.(|TT|), the thermal fluctuation force, n = term, 

2dqs 



in this case is 
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[(2dg,)4e2<i,.^^(2dq,)] . (16) 
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where Ei{x) is the exponential integral. Again as g'gd ^ 
1, the force between atom and graphene becomes /t = 
3/cBTx(0)/(4(i^), which is the thermal interaction between 
an atom and a perfectly conducting plate. 

The normalized forces, F/Ft and ///t, as a function of 
2dqs are shown in Fig.([2|. One finds that for large dis- 
tances spatial dispersion becomes less important, and the 
force goes to the one expected for ideal metallic surfaces. 
For short distances the electron screening shields the force 
such that it goes to zero at very small distances. From 
Fig.( |2(a)"] ) it is seen that dispersion is significant when 
^5^-10 or less. For the case of gapless graphene, this 
corresponds to a distance of d ^ 20 nm at room tempera- 
ture, in which case higher Matsubara frequencies need to 
be considered. Hence, for the case of graphene, dispersion 
does not play a significant role in the region where the 
fluctuation force is thermal. Spatial dispersion becomes 
important for graphene at larger distances, when qs de- 
creases. 

It has been noted that while for good conductors and 
dielectrics thermal Van- der- Waal /Casimir forces are not 
heavily influenced by dispersion effects, the latter may 
have a strong influence on the force for poor conductors 
[31]|32]. A similar result occurs for graphene, where upon 
the conduction properties of the material are decreased by 
introducing a band gap. 

By inducing an energy band gap, the effect of spatial 
dispersion in graphene is increased by effectively decreas- 
ing ^s, as shown in Fig.( |2(b)] ). It is clear that the influence 
of dispersion increases as qs gets small with the effect that 
the thermal fluctuation forces between graphene sheets are 
substancially decreased. For example for A = 0.05 eV ^ a 
typical value for the gap, one obtains a distance for disper- 
sion to have strong effects of d 40 nm, within the range 
where thermal fluctuation forces are dominant. Of course, 
by inducing a sufficiently large band gap, thermal effects 
cease from becoming dominant. This effect is analogous 
to 3-D materials, whereupon the influence of dispersion in 
the fluctuation forces of a conductor become strong were 
it to become a poorer conductor. 
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Fig. 2: (a) Thermal stress between two graphenes normalized 
to the thermal stress between two metals and thermal force 
between graphene and atom normalized to the thermal force 
between metal and an atom, (b) qs normalized to the value at 
zero band gap ^so as a function of the energy band gap A for 
different temperatures, T ^ 2m K , T ^ ?>m K , T ^ Am K . 



T = K. When the temperature is zero, the sum- 
mation in F (Eq.([8|) is substituted by an integral, = 
[h/{27ikBT)]j^ duj. The T = K case corresponds to 
purely quantum mechanical contributions to the Casimir 
stress. If the graphene conductivity is given by the uni- 
versal value (Jo = e^/(4fi), the graphene/graphene Casimir 
force per unit area is found as ^15|jT7] 



(17) 



where Fq = ii'^hc/ {24Qd^) is the force per unit area be- 
tween two perfectly conducting metallic sheets and i<i{a) = 
720a/(327r^) with a = 1/137 being the fine structure 
constant. It is interesting that the force in this limit 
does not depend on any quantum mechanical character- 
istics nor the speed of light. Comparing Fg to Fq shows 
that the distance dependence of the graphene/graphene 
force is the same as for metallic plates in the retarded 
regime, however, the magnitude is reduced by an order of 
the fine structure constant. Incidentally, this differs from 
the non-retarded "2D" metal and doped graphene result- 
F - (^-^/2 [MlfT? 



The presence of /i or A alters the graphene conductiv- 
ity, which can change the Casimir interaction significantly. 
Using the isotropic cro{iuj) from Eqs.([6j[7|), it is found that 
F/Fq increases almost linearly as a function of d when 
/i 7^ 0, while for A 7^ and /i 7^ the normalized Casimir 
stress has nonlinear dependences - Fig. (3(a)). This hap- 



pens because at shorter distances the effect of the band 
gap is more pronounced, while at larger distances the fi- 
nite chemical potential has a stronger effect. It is found 
that F/Fq can be diminished by increasing the band gap 
of the graphene sheets ( Fig.( |3(b)] )) or it can be enhanced 
by increasing the chemical potential ( Fig. ( 3(c)'] )). These 
changes can be quite substantial, providing that the sys- 
tem supports large /i or A. 

It is interesting to consider to what extent spatial dis- 
persion affects the graphene/graphene Casimir interac- 
tion. In the absence of a band gap or chemical poten- 
tial, dimensional analysis dictates that the force may be 
written as 



F = Fog{a,vo/c), 



(18) 



where the scaling factor g{a,vo/c) now depends on the 
fine structure constant, and the ratio of vq/c. If one uses 
Eq.(|4|, the scaling factor is evaluated as f{a,vo/c) ~ 
0.0043. Note that the distance dependence is still 
F ~ Comparing the scaling factors tz and g 

shows that the inclusion of spatial dispersion reduces the 
graphene/graphene Casimir attraction by ~ 20 %. If 
one were to approximate the conductivity of graphene 
via RPA, then inclusion of dispersion would reduce the 
graphene/graphene interaction by ^ 9 %. Figs.( |3(b)|3(c)| ) 
further show that dispersion effects are rather insignificant 
as /i and A increase. 

T 7^ OK. Finally, the full finite temperature Casimir 
interaction is considered. Eq.([8| indicates that the largest 
contribution to the force occurs when qd ^ 1. On the other 
hand it was shown that for finite frequencies, the conduc- 
tivity varies sharply as a function of wave number when 
qs > 10^ cm~^ - Fig(l(b)). Therefore, dispersion effects 
start becoming important for Matsubara terms in Eq.([8| 
greater than zero at a distance of ds ~ l/qs = 0.01 pm. 
Fig. (4(a)) shows the Casimir interaction force per unit 
area between two graphene sheets where spatial disper- 
sion is only included in the force for the n = thermal 
term and also where spatial dispersion is included for all 
Matsubara frequency terms. If the distance is reduced 
much further the force reduces back to the zero tempera- 
ture result. 

Much more influential parameters to the Casimir inter- 
action are the induced chemical potential or energy band 
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Fig. 3: (a) Casimir stress normaUzed to Fq at T=0 and cal- 
culated using the isotropic conductivity ao{i(jj) for different A 
and /i. The flat line indicates the normalized stress obtained 
via the universal conductivity ctq; (b) Casimir stress normal- 
ized to Fq at T=0 as a function of A calculated using ao{i(jj) 
and a{iuj,q) (/x = 0); (c) Casimir stress normalized to Fq as a 
function of /a calculated using ao{i(jj) and a{i(jO,q) (A = 0) 



gap. The effects are displayed in Fig.( |4(b)] ), which shows 
the stress between graphenes for different band gaps, and 
chemical potentials. In addition, the role of the chemical 
potential is to increase the magnitude of the force while 
the energy band gap decreases it (not shown graphically). 
The characteristic behavior is similar as for the quantum 
T = K limit. 



Fig. 4: (a) Casimir stress normalized to the idealized stress 
Fo at T = 300 K with and without the inclusion of disper- 
sion for Matsubara terms greater than zero, (b) Casimir stress 
normalized to the idealized stress Fq Sit T — 300 K. 



conductivity. Spatial dispersion corrections start appear- 
ing at distances of d ^ 0. 01 /im at room temperature. 
A finite chemical potential can enhance the Casimir in- 
teraction substantially and further inhibit spatial disper- 
sion. If a band gap is induced at finite temperature, how- 
ever, graphene exhibits similar characetristics to a poor 
metal and dispersion effects start becoming more impor- 
tant, that is the thermal fluctuation forces are reduced by 
the presence of electronic charge screening. 

In the present work, dispersion has been investigated 
using the 2-band model. Yet at shorter distances, in the 
order of 10 nm or less, higher band transitions become im- 
portant for the Van-der Waals/Casimir interaction. The 
inclusion of such effects without spatial dispersion have 



already been investigated in a previous work 10 
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Conclusion. — The Casimir interaction between two 
graphene sheets has been investigated at room and zero 
temperatures by using linear response theory with a con- 
ductivity response function in order to determine how spa- 
tial dispersion affects the interaction. It was found that 
these effects are not of major importance, either at large 
or low temperatures. The Casimir interaction may be 
well approximated by using a frequency only dependent 
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